Expansion of a Finite Size Plasma in Vacuum 
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The expansion dynamics of a finite size plasma is examined from an analytical perspective. Results 
regarding the charge distribution as well as the electrostatic potential are presented. The acceleration 
of the ions and the associated cooling of the electrons that takes place during the plasma expansion 
is described. An extensive analysis of the transition between the semi infinite and the finite size 
plasma behaviour is carried out. Finally, a test of the analytical results, performed through numerical 
simulations, is presented. 
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The production of energetic particles through the in- 
teraction of an intense laser pulse with a solid target is 
a topic that has been widely investigated during the last 
three decades. For sufficiently long pulses, the emitted 
particles emerge from a coronal plasma formed by the 
laser on the target foil and the expansion of the plasma 
in vacuum plays a key role in such a phenomenon. Thus, 
despite in Refs.Q, 0- 03- EH the expansion of a semi infinite 
plasma has been widely investigated, present experiments 
often involve "thin foils" as targets [1 S SS E E3, E3 
which are in some cases no thicker than a few tens of 
a Dcbyc length, and we believe that, in these cases, the 
experimental results must be analyzed in terms of the 
expansion of a finite size plasma. The finite-size and the 
semi-infinite cases are very different from each other, as in 
the latter case an infinite amount of energy is available. 
As a consequence, even if the ions are accelerated, the 
energy of the electron plasma remains constant. On the 
contrary, in the case of a finite size plasma an exchange 
of kinetic energy between electrons and ions takes place 
as long as the two populations have different velocity dis- 
tributions, i.e., as long as a charge separation is present. 
The aim of this work is to provide a detailed analytical 
description of the thermal expansion of a globally neutral, 
finite size, unidimensional plasma. Analytical predictions 
will be compared with numerical results obtained with a 
Particle in Cell (PIC) code. 

The configuration at the initial time to, which formally 
can be defined as the hydrostatic equilibrium in the limit 
of infinitely massive ions, is specified by the ion den- 
sity riio(x) = ni(x,to) = uq 9{a — |x|), and by the 
Boltzmann-like electron profile n e o{x) = n e (x,to) = 
n exp (e$(x, to)/T e o), where a is the half-thickness of the 
plasma, 2ano is the total number of positively charged 
particle, 8(x) = for x < and 6(x) = 1 for x > 0, 
$(x,to) is the electrostatic potential, n is the density 
at the position where &(x,to) = 0, £e0 the initial elec- 
tron plasma temperature and n e o(x)dx = 2aiiQ be- 
cause the system is globally neutral. For the sake of 



notational simplicity we are considering a plasma where 
the ion charge is equal and opposite of that of the elec- 
trons. Measuring space in units of the initial electron 
Debye length A^.o = (T e o/47re 2 no) 1 / 2 and rescaling the 
potential <&(x,t) as $(x,t) — > 4>{x,t) = e$(x, t)/T e0 , we 
rewrite the Boltzmann equation in the simpler form 



n e o(C) 



ne 



HO 
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In order to obtain the electrostatic potential $(x, to), 
and hence n e o(x), we need to solve the Poisson equation 
V 2 $(x,t ) = — 4:Tre[n i0 (x) — n e0 (x)], coupled to Eq.Q. 
Because of the inversion symmetry (£ — > — £) of the con- 
figuration, we will restrict our analysis to the half-plane 
C > 0. We choose the potential to be zero in £ = a, where 
a denotes now the normalized half width of the plasma 
slab, so that n = n e g(a) = n a and Poisson equation reads 



=0(0 



no 



1 6(a-0+( ^e*<0) 6((-a) 



no 



(2) 

The different functional form of riio(x) for |ar| < a (inter- 
nal region) and for \x\ > a (external region), implies that 
separate treatments are required. First we will calculate 
the electron density at the ion front x = a. Using this 
result, we will obtain an approximate analytical expres- 
sion for the electrostatic potential in the internal region, 
while in the external region we will use the analytical so- 
lution of Eqs.(2[2J found in @. 

The value of n a can be derived in the following way. In- 
tegrating Eq.J2J once gives the electric field E(Q in the 
form 

£7(0 = y/2[(n a /no) (e*«) - e*«») - (</>(C) - 0(0))] (3) 
for < C < a and, for £ > a, in the form 



£7(C) = ^2(n a /n )e*(C), 



(4) 
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where the integration constants have been fixed such that 
in Eq.<J3|) £7(0) = 0, as follows from the inversion symme- 
try, and that in Eq.Q lim^oo £7(0 = which expresses 
charge neutrality (together with lim^_ too (n a /no)e*^^ = 



Rarefaction front 1Z x f(t) 



FIG. 2: Model profile of the ion density as plasma expands. 



FIG. 1: Behaviour of n a /no and of n e o(0)/no as a function 
of a. 



since there are no electrons at infinity). Imposing the 
continuity of E(£) across the boundary £ = a and us- 
ing Eq.QJ, we obtain the following relationship for the 
potential and the electron density at £ = 0(0) = 
(n a /n )e* (0) = n e0 (0)/n . The integration in Eq.Q for 
</>(x) can not be performed explicitly. Numerical integra- 
tion gives the electron density n e o(0) at ( = and the 
electron density n a at £ = a as a function of the plasma 
length a. Fig. shows that, for a > 1, (i.e., in dimen- 
sional units for a ^> Ad,o)j we find (f)(0) = n e o(0)/no ~ 1 
and n a /no ~ 1/c, so that for most cases of interest, the 
electron density at the ion front is equal to that obtained 
in the limit of a semi-infinite plasma. The solution of 
Eq.Qj with the condition 0(a) = 0, is given by 



E = 



2\/n a /n 



V 'n a / 'n (C - a) + a/2' 
with electrostatic potential given by 

<t>(() = In 



y/n a /n (C — a) + n/2^ 



(5) 



(6) 



In order to derive an approximate solution of Eq.©, we 
define a layer of thickness (5 inside the plasma, specified 
by the assumptions that n e o(C) — no for < C < (a — 
where S ~ and ~ const = 0(0), and that the electron 
density distribution differs significantly from the fixed ion 
distribution only for £ > (a — ft). Inside this layer we 
adopt an approximate parabolic fit for the electrostatic 
potential of the form 0(C) = d(( - a) + C 2 (C - a) 2 A 
where the constant term has been set equal to zero since, 
by continuity with the external solution, 4>(a) = 0. This 
fit of the potential corresponds to a rough approxima- 
tion in Poisson's equation where the electron den- 
sity inside the layer is taken to be constant and equal 
to an intermediate value between its two values, n a and 
no, at the borders of the layer. This value must be de- 
termined, together with G\ and 8, sclf-consistently by 
requiring that the potential and the electric field be con- 
tinuous at £ = a, a — 5, which ensures in particular charge 



conservation. From the continuity of the e lectric field at 
£ = a we obtain from Eq.JSJ) C\ = —y/2/c. At ( = a — 5, 
the continuity of the electric field gives C2 = —C\/6, while 
the continuity of 4> defines 5 as S = s/2e(/>o ps s/2e. 
We now derive a physical model capable of describing 
the electron cooling that occurs as the plasma expansion 
takes place and, consequently, ion acceleration. Measur- 
ing temperatures in units of T , where T = T e0 , mass in 
units of the ion mass m^, and introducing the dimension- 
less time unit r = w p it, with lo v \ = \/4Tmoe 2 , we denote 
with T e the electron time dependent temperature, with 
m e the electron mass and assume that both the initial 
fluid and thermal ion energy is zero. The quantity T e 
can in principle depend both on x and r, however in 
what follows it will be assumed spatially constant. In 
order to derive an analytical expression for the cooling 
of the electron population we need to introduce a sim- 
plified description of the ion expansion and of the time 
evolution of the electric field. For this purpose, we in- 
troduce the approximate ion profile nj(x, r) sketched in 
Fig-© which is specified once the quantities x s and xj 
are known. The point x s models the position of the rar- 
efaction front, which obeys the equation of motion 



= a — 5 — 



s dr. 



(7) 



with c s the time dependent ion acoustic speed which, in 
dimensionless units, is given by c s = \/T~ e . The quantity 
Xf represents the ion front position, defined by the con- 
dition rii(x, r) = for x > Xf, which obeys the equation 
of motion 



'if = E f, 



(8) 



with initial conditions i;/| T= o = 0, Xf\ T= o = a and Ef 
the self consistent electric field at point x = Xf . Profile 
then becomes nj(x,r) = n 8(x s — \x\) +n D 9(\x\ — 
x s ) — noD 6(\x\ — Xf), with D = (a — x s )/(xf—x s ), which 
corresponds to redistributing the charge no (a — x s ) on 
the interval (xf — x s ) 01 . Assuming hydrostatic equi- 
librium at every time r between the electron and ion 
plasma we find D n^/nf = c for the electron den- 
sity rif at point x = Xf, thus Ef = \J (2/e) DT e , which 
completes the description of both the ion profile and dy- 
namics. Regarding the electrons, we assume quasineu- 
trality in the whole region [0, Xf — S] except for a small 



3 



interval [ 1, with Ax to be determined, in 

which charge separation occurs such that Boltzmann re- 
lation <£Q) can be satisfied. The electric field E cs in 
[x s — Ax, x s ] is approximated by E cs ~ (1 — D), which is 
the field of a plane capacitor with surface charge den- 
sity cr = (1 — D) /At: in dimcnsionless units. There- 
fore - f*°_ Ax E cs dx = T e [(j) (x s ) -<f>(x s - Ax)} which 
defines Ax as Ax = -T e (lnD)/(l - D). The cooling 
phenomenon is associated with the reflections that the 
electrons undergo at the potential barriers in the re- 
gion [x s — Aa;,3; s ] (internal barrier) and [xf — S, +oo] 
(external barrier), which are moving with velocities — c s 
and Vf, respectively. Assuming the reflections to be 
elastic, m e <C 1 implies that the power loss due to 
a reflection of an electron with velocity v is given by 
AU(v) = 2m e (v — v)v, with v = {— c s , Vf}. The corre- 
sponding power is given by AU(v) times the number of 
electrons with velocity v > v that hit the barrier per unit 
time, integrated over the electron distribution function. 
Assuming the latter to be Maxwcllian, the number of 
electrons with velocity in the interval [v, v + dv] reflected 
per unit time by the internal and external barrier is given 
by (1 — D)novf(v) and Dn$vf(v) 9(v — Vf) respectively, 
where f{v) = (l/V^ v th ) e~ v '' V 2 **" and v th = y/T e /m e . 
The total energy of the electron plasma is U = noaT e /2, 
hence, taking the temperature T e to be uniform, dU/dr = 
(nQa/2) dT e /dr and, taking Vf, c s <C v t h, we find 

^r-^- = -(Vo.)[Dv f -(l-D)c a ], (9) 

which is the electron cooling equation, to be coupled 
both to Eq.Q and Eq.(J5J) and then integrated numer- 
ically. The time dependent total kinetic energy Ui of 
the ions satisfies the conservation equation Ui = noa/2 — 
noaT e /2 — U e i, where U e i is the contribution from the 
release of the electrostatic energy present in the initial 
"equilibrium" configuration described previously. Note 
that 2U e i/{T e ) w (2dU e i/dt) / (dT e /dt) « S/a and thus 
Ui ~ 1 — T e /2 except for the very beginning of the expan- 
sion process, when the electrostatic energy contribution 
is significant and must be taken into account. 
We may summarize the results of the above modeling 
of the plasma expansion in vacuum by referring to a, 
r/(a^/m e ) and r/y/a as the relevant physical quanti- 
ties characterizing the plasma dynamics. The condition 
a » 1 implies that, at every time t, the electron den- 
sity at the expanding front is rif = Dn^/e, as in the 
hydrodynamic equilibrium of a semi-infinite slab. The 
quantity r j (a-Jm e ), where a^/m e is the typical time it 
takes an electron to cross the target, roughly gives the 
ratio of particles which, at time r, have undergone at 
least one reflection at the ion front. Thus the initial 
expansion phase when r j (a,y/m e ) <C 1 corresponds to 
the semi-infinite, constant temperature plasma limit dis- 
cussed in dHla. The quantity t/^E measures, at the 
early stage of the process, the expansion of the ion com- 
ponent in the plasma. In fact, taking Vf oc r, which 
corresponds to a uniformly accelerated motion of the 
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FIG. 3: Plasma expansion regimes in the (a, r)-plane. Elec- 
tron cooling is important in region C. The uniform electron 
temperature assumption used in our analytical model applies 
for a < l/m e . 




FIG. 4: Total normalized ion kinetic energy versus time r as 
obtained from the analytical model (dashed line) and from a 
PIC simulation (solid line) in a plasma with a — 50 A^. 



front and neglecting coefficients of order unity, we obtain 
(xf — a) /a = dV/V as r 2 /a where V is the volume oc- 
cupied by the ions. In the initial expansion phase, from 
Eq.© we obtain dT e /T e as -2 (v f T)/a as -2 dV/V, 
which implies that in this phase the plasma behaves like 
a perfect monoatomic gas obeying the law of unidimcn- 
sional adiabatic expansion. Thus, the plane (a, t) can 
be divided into four regions by the two curves r = ^/a 
and r = a^/m e , as shown in Fig.Q. The two curves 
cross at a = l/m e . In the three regions A,B,D defined 
by r < max(ay/m e , y/a), electron cooling is not impor- 
tant and the slab expansion can be approximated by the 
semi infinite case. In region C the effects of the the finite 
size of the plasma becomes apparent and ion acceleration 
takes place together with electron cooling. 
Now we reexamine the approximations used in the deriva- 
tion of Eq.@. The assumption of uniform temperature 
rests on the comparison between the energy redistribu- 
tion time, approximated by a^Jm e , and the typical ion 
plasma expansion time y/a, which roughly represents the 
energy confinement time. Therefore the temperature in- 
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FIG. 5: Normalized electric field Ej at the ion front (a) and 
ion front position xj (b) versus time r as obtained from the 
analytical model (dashed line) and from the PIC simulation 
(solid line), in a plasma with a — 50 \d- The dotted line 
corresponds to semi-infinite case (J. 



side the plasma can be taken to be uniform as long as 
am e <C 1 (the area to the left of the vertical clotted line 
in Fig.(j2Jl) and the predictions of our analytical model 
are thus restricted to this domain which, however, covers 
the typical experimental configurations H, [n| . The ap- 
proximation that the electron reflections are elastic in the 
frame co-moving with the reflecting barrier, can be writ- 
ten as y/m e d[\n (T e <f>)]/dT <C 1, where we take the elec- 



tron reflection time to be of order ^/m e . The validity of 
this condition has been verified by integrating Eqs.QIHl 
E5J) numerically. Finally, because of the significant spread 
of the ion front, the effective value of Ef is smaller than 
the one obtained within this model. This may be cor- 
rected by introducing the reduction factor 1 — (1 — D) a , 
with a a constant to be fixed numerically. In order to 
determinate the value of a and to verify our analytical 
results, a series of numerical simulations has been per- 
formed using a unidimensional PIC code. The numerical 
results that we present here describe the expansion of a 
plasma slab 50 X c i thick with initial conditions given by 
Eas. HH6|) . An overall good agreement between the ana- 
lytical and the numerical results is achieved with a = 1.4. 
The ion total kinetic energy and the value of Ef and 
the position of the ion front are plotted versus time in 
Figs. and JSJ), respectively. Different simulations for 
slabs up to 150A,i thick have been performed and have 
confirmed the validity of the analytical model and in par- 
ticular of the choice of the value of the parameter a. The 
comparison between the finite size results and those de- 
scribed in pj is shown in Fig.© and indicates that the 
semi-infinite and the finite size plasma behaviours are 
very similar only at the early stage of the expansion pro- 
cess, while they become, as expected, substantially dif- 
ferent when the electron kinetic energy is significantly 
reduced. 

In conclusion our analytical model shows that, in the case 
of a finite size plasma, the process of electron cooling 
leads to a nearly complete energy transfer from electrons 
to ions on time intervals of the order of a few tens of the 
transit time a^/me. This sets an upper limit to the en- 
ergy ions can acquire in the thermal expansion of a finite 
size plasma. Furthermore, different expansion regimes 
as function of time and plasma size have been identified 
and the numerical verification of the proposed model has 
been discussed. 
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